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Abstract 

Previous works have shown numericahy that the response of a "stochastic res- 
onator" is enhanced as a consequence of spatial couphng. Also, similar results 
have been obtained in a reaction-diffusion model by studying the phenomenon 
of stochastic resonance (SR) in spatially extended systems using nonequilib- 
rium potential (NEP) techniques. The knowledge of the NEP for such sys- 
tems allows us to determine the probability for the decay of the metastable 
extended states, and approximate expressions for the correlation function and 
the signal-to-noise ratio (SNR). Here, exploiting known forms of the NEP, we 
have investigated the role of NEP's symmetry on SR, the enhancement of the 
SNR due to a selectivity of the coupling or diffusion parameter, and discussed 
competition between local and nonlocal (excitatory) coupling. 
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I. INTRODUCTION 



Many recent theoretical and experimental studies have provided a large body of evidence 
on the essentially constructive role played by fluctuations in a variety of intriguing noise in- 
duced phenomena. Some key-examples are: problems of self-organization and dissipative 
structures, noise induced transitions, noise induced phase transitions, thermal ratchets or 
Brownian motors, coupled Brownian motors, noise sustained patterns and stochastic reso- 
nance 

This last phenomenon, that is stochastic resonance (SR), has attracted considerable 
interest due, among other aspects, to its potential technological applications for optimizing 
the transmission of information in nonlinear dynamical systems. For a comprehensive recent 
review see Ref. [^, that shows the large number of applications in science and technology, 
ranging from paleoclimatology, to electronic circuits, lasers, and noise-induced information 
flow in sensory neurons in living systems, to name a few. Several recent papers have aimed 
at achieving an enhancement of the output SNR by means of the coupling of several SR 
units in what conforms an "extended medium" 

The relevance of pattern-formation phenomena to several areas of science and technology 
is a very well established fact. Accounts of the state of the art can be found in many reviews 
and books collecting the results obtained over the last couple of decades regarding the 
description of pattern formation and propagation phenomena in self-organizing systems. 

It is a common belief that the nontrivial spatio-temporal behaviour occurring for instance 
in the Complex Ginzburg-Landau Equation (CGLE), reaction-diffusion (RD) schemes, and 
other nonequilibrium systems, originates from the non-potential (or non-variational) char- 
acter of the dynamics, meaning that there is no Lyapunov functional (LF) for the dynamics. 
However, Graham and co-workers, have shown that a Lyapunov-like functional does exist 
for the CGLE 0. 

In order to flx ideas, we start by discussing the existance of LF in two different dynamical 
situations. The simplest case, in which a LF exists, corresponds to a gradient flow system 
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such as 

d 

Xj = - —V{xi, ..,x„). (1) 

The fixed points will correspond to the extrema of the "potential function" V{xi, ..,Xn), 
and the system will evolve towards these minima of V{xi, following trajectories cor- 

responding to the line of steepest descent. Clearly V{xi, .. LF as it also satisfies 

dv ^ dvdx^ ^ -vf^V < n (2) 

dt ^ dxj dt ^ [dxjj - ' ^ ^ 

Now we consider an example corresponding to a non-gradient fiow. We take the following 

case 

dV 



X 



where T is an arbitrary, positive definite matrix. We can separate it into a symmetric {S) 
and an antisymmetric {A) part 

T = S + A 
S = -{T + T^), S =S^ 

A=\{r - T^), A = - A^. (4) 

The fixed points are given by the extrema of V . On the other hand we have that V also 
fulfills 

dV dV dV dV dV 

H = - Y.(-^hg^^g^^ < 0. (5) 

as, clearly, the first term on the rhs is < 0, while the second one is = 0, hence, y is a LF. 
The system evolves to the minima following trajectories different from the steepest descent 
ones, determined by S, since the antisymmetric part of T induces a flow in the system that 
keeps the LF constant. A thorough (and didactic) discussion on the LF and different kinds 
of dynamical flows can be found in 0. 

The use of the LF concept (as well as the related one of nonequilibrium potential for 
stochastic systems) in relation to far from equilibrium systems is scarce. However, there are 



recent papers that have made use of these ideas in relation not only with CGLE but 
also for reaction-diffusion systems. Some of these works concerned the derivation of effective 



equations for the evolution of fronts ||rT|; while others were related to the determination of 
the global stability of the resulting patterns and the possibility of exchanging the relative 
stability between attractors ■ 



In the next Section we introduce the concept of nonequilibrium potential (NEP). When 
the NEP can be obtained, such an approach offers an alternative way of confronting a 
problem that has attracted considerable attention, both experimentally and theoretically. 
Namely, the relative stability of the different attractors, corresponding to spatially extended 
states, and the possibility of transitions among them due to the effect of fluctuations . The 
latter aspect which is of great relevance in the study of SR in spatially extended systems, is 
the objective of this work. 

The organization of the paper is as follows. As indicated, in the next section we briefly 
discuss some basic notions about nonequilibrium potentials, and show a few relevant exam- 
ples for reaction-diffusion systems. In Section III we present the results for the SNR in some 
of the previous examples. Section IV contains a final discussion and some conclusions. 



II. NONEQUILIBRIUM POTENTIAL 
A. Brief Review 

Loosely speaking, the notion of non-equilibrium potential (NEP) corresponds to an exten- 
sion of the notion of equilibrium thermodynamical potential to non-equilibrium situations. 
In order to introduce such NEP, we consider a general form of non-linear stochastic equa- 
tions, admitting the possibility of multiplicative noises. In particular, we consider equations 
of the form 

r = K''{q)+g'^{q)Ut), z/ = l,...,n; (6) 
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where repeated indices are summed over. Equation (|^) is stated in the sense of Ito. The 
{i = l,...,""^ < n)} are mutually independent sources of Gaussian white noise with 
typical strength rj. The Fokker-Planck equation corresponding to Eq.(^) takes the form 



dP d n d'^ 

d,t dq'^ 2dq''dq''^ 



(7) 



where P{q,t;ri) is the probability density of observing q = {qi,---,qn) at time t for noise 
intensity rj, and Q'^^{q) = di (q) di (q) is the matrix of transport coefficients of the system, 
which is symmetric and non-negative. In the long time limit {t oo), the solution of Eq.(^ 
tends to the stationary distribution Pstat{(l)- According to [P], the NEP associated to 

Eq.(^, is defined by 



= -limr] In Pstat{q,v)- 



(8) 



In other words 



Pstat{q)d'^q = Z{q) exp 



7] 



+ 0(r/) 



9' 



where is the NEP of the system and Z{q) is defined as the limit 



InZ(g) = lim 



In Pstat{q,v) + -^{q) 

7] 



Here dQ,, 



is the invariant volume element in the g-space and G{q) is the determinant 



of the contravariant metric tensor (for the Euclidean metric it is G = 1). It was shown 
that $(g) is the solution of a Hamilton- Jacobi like equation (HJE) 



9$ 1 9$ 9$ 



(9) 



and Z{q) is the solution of a linear first-order partial differential equation depending on $(g) 
(not shown here). 

Equation (^) and the normalizability condition ensure that $ is bounded from below. 
Furthermore, from the separation of the streaming velocity of the probability flow in the 
steady state into conservative and dissipative parts, it follows that 



^ = ^ ^^^^ = "2 ^ 9^ 9^ - °' 

i.e. $ is a LF for the dynamics of the system when fluctuations are neglected. Under the 
deterministic dynamics: = K'^{q), $ decreases monotonically and takes a minimum value 
on attractors. In particular, $ must be constant on all extended attractors (such as limit 
cycles or strange attractors) 

B. Examples of Nonequilibrium Potentials 

1. Scalar Reaction-Diffusion Model 

As a first example we focus on a one-dimensional, one-component model of an elec- 
trothermal instability [|l|, which corresponds to an approximation to the continuous limit 
of the coupled system studied by Lindner ei a/ [Q. For this model, the effect of boundary 
conditions (b.c.) in pattern selection, the global stability of nonhomogeneous structures, and 
the critical-like behaviour due to the coalescence of two patterns [§, have been studied. 

The RD model that we work with describes the time evolution of a field t) which 
represents the temperature profile in the so-called "hot spot model" in superconducting 
microbridges (or ballast resistor) The evolution of (p is given by 

^ = D - + ^[0 - 0J + e(a;, t), (10) 

where ^(x, t) is a white noise in space and time, in the bounded domain x G [—L,L] with 
Dirichlet b.c. at both ends, i.e. (f){±L,t) = 0. 6[x] is the step function. We restrict our 
discussion to the parameter range where the (associated deterministic) system has two stable 
attractors (patterns). The piecewise linear approximation of the reaction term, mimicking 
a cubic polynomial, was chosen in order to find analytical expressions for the spatially 
symmetric solutions of Eg . ([10|) . It is clear that the trivial solution (j)o{x) = 0, which is 
linearly stable, exists for the whole range of parameters. Besides this solution we find only 
one stable nonhomogeneous structure, 0s (a;), which presents an excited {4>s{x) > 0c) central 
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zone, and another similar unstable structure, with a smaller excited central zone. The 

latter pattern corresponds to the saddle separating both attractors (f)o{x) and (psi^) P,p^|. 



The indicated patterns are minima of the NEP of our system. For the present case, such 



a NEP reads ||T0| 



It can be shown that |^ = — |^ and also $ = — / (^^^ dx < 0. This functional offers us the 
possibility to study both the local and global stability of the patterns as well as the changes 
associated to variations of model parameters. p|,p!0| 

In Fig. 1 we depict $[0,0c] evaluated at the stationary patterns 0o ('^'[0o] = 0), 0s(x) 
(^<|)* = $[0^]) and (pui^) = for a system size L = 1, as a function of (pc for a 

given value of D. The upper branch is the NEP for where $ attains an extremum 

(as a matter of fact it is a saddle). On the lower branch, for 4>s{x), and also for 4>o{x), the 
NEP has local minima. The curves exist up to a certain critical value of 0c at which both 
branches collapse P,p!0|. It is interesting to note that, since the NEP for 0m (x) is always 
positive and, for 0s (x), > for "large" values of 0c and also < for "small" values 
of 0c, vanishes for an intermediate value 0c = 0*, where 0s (x) and 0o(x) exchange their 
relative stability. 

2. Three Component Activator-Inhibitor Model 

Here we present an exact form of the nonequilibrium potential for a three component 
system of the activator-inhibitor type (with one activator and two inhibitors) in a particular 
parameter region. Such a three component system provides the adequate framework for 
a minimal model describing pattern formation in high-pressure or low-pressure chemical 
reactors [0. The system we consider is described by the following set of equations 



^^^^ = DVM^, t) + f{u{x, t)) - v{x, t) - w{x, t) + g\Ux, t) + gl^U^, t) 
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^^dv{x^ = u^VM^, t) + (3u{x, t) - jv{x, t) + g'l^.ix, t) + g^U^, t) 



£2 



dt 
dw{x, t) 

m 



U2'^^wix, t) + P'u{x, t) - ywix, t) + (7rei(a;, t) + g^Ux. t), 



(12) 



where x represents an n-dimensional spatial coordinate. The t)'s are gaussian white- 
noise sources of zero mean satisfying 



(13) 



where r] is again a small parameter measuring the noise intensity. All the parameters and 
fields shall be considered as dimensionless (scaled) quantities. We shall only consider situa- 
tions where the noise terms in the third of Eqs. (0) are negligible, and we set gi = g2 = 0. 
Furthermore, we will analyze the system in the limit ui = €2 = with ei = 1 and 1^2 = 
where for the now temporally slaved inhibitor w we have 



w{x,t) = dx'G{x, x')u{x' ,t), 



(14) 



where G{x,x') is the Green function of the third of Eqs. (]T^) in the indicated limit [§. 
Hence the system can be reduced to an effective two component system with a nonlocal 



interaction 13.12 



When the relation 
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Qu(3 + Q^ 



(15) 



where Q„ = «)^ + (^2)^ Qv = (^i)^ + (^2)^ and = ^"^1 +^2^2^ holds, the two equations 
for such an effective two component system adopt the form 



/ 


du{x,t) 


] 






( 




at 




= -T 


Su(x,t) 


+ 




dv{x,t) 


) 




5<}> 




V 


dt 




\ 5v{x,t) / 


V 



glUx,t) + gl^Ux,t) 



(16) 



Here the matrix T, which has the form of the matrix T in Eq. (|^) is given by 
1 



r 





_ 1 


(q q \ 


+ 


n Quv ^ 
^ 2 


Qv J 


~ 2 


^ Quv Qv J 




Quv n 



S + A, 



(17) 



and the functional ^[u{x,t),v{x,t)] by 
^[u{x,t),v{x,t)] = dx 



D If 

— (^u{x,t)Y + V{u{x,t)^v{x,t)) + -— I dx'G{x,x')u{x,t)u{x' ,t) 

(18) 



where 



y{u,v) = -— f[u)du + ——u + —V -2—uv. (19) 

When the symmetric matrix S is positive definite (when QuQv > Qtv holds), the func- 
tional $, that in the associated deterministic system decreases monotonously with time, is 
the NEP of the system in Eqs. (|TBp satisfying the HJE (Eq. (|^)). Equation (|TB|), which 
resembles a detailed balance condition, arises in this context as a mathematical constraint 



necessary for $, as defined by Eq. ( |T8|) , to be the solution of the HJE above mentioned. 

We limited the analysis to the parameter region where Eq. (^) is valid, the matrix S is 
positive definite, and hence the nonequilibrium potential is given by Eq. ([T8|) . Although these 
conditions impose restrictions on the range of application of our treatment, it is worth noting 
that, after choosing the values of the (7j"s satisfying the condition of positive definiteness. 



there is still a wide spectrum of interesting situations to analyze [12,|T3||. Furthermore the 



nonlinear function f{u) is still arbitrary, opening a wealth of possibilities. When plotting 
$vs.0c) with $ evaluated on the stationary patterns, we see a result similar to the one shown 
in Fig. 1. 

III. STOCHASTIC RESONANCE IN EXTENDED MEDIA 

A. Preliminaries 

The calculation of the SNR proceeds, for the spatially extended problem, through the 
evaluation of the space-time correlation function {(f){y,t)(j){y' ,t')) . To do that we have used 
a simplified point of view, based on the two-state approach |T^, which allows us to apply 



some known results almost directly. We consider a random system described by a discrete 
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dynamical variable x adopting two possible values: ci and C2, with probabilities ni^2{t) re- 
spectively. Such probabilities satisfy the condition ni(t)+n2(t) = 1. The equation governing 
the evolution of ni{t) (with a similar one for n2{t) = 1 — ni{t)) is 

drii dn2 



dt dt 



W2{t)n2{t) - Wr{t)n,{t) = W2{t) - [W2{t) + W,{t)]n,, (20) 



where the Wi^2{i) are the transition rates out o/the x = C\^2 states. For the bistable extended 
system in which we are interested, such states correspond to the spatially extended attractors 
(see Refs. p|,p!6| for details). 



If the system is subject (through one of its parameters) to a time dependent signal of 
the form Acos(co'st), up to first order in the amplitude (assumed to be small) the transition 
rates may be expanded as 

W\{f) = /il — ttiA cos(co'st) 

W2 (t) = /i2 + 02^ COs(cJ,t) , (21) 

where the constants /xi,2 and a\f2 depend on the detailed structure of the system under study. 
Here we remark that the /ij's, which are the (time independent) values of the PVj's without 
signal, are in general different from each other as a consequence of the different stability 
of the two states, and the same happens to the a^'s. With the indicated modulation the 
system becomes nonstationary but we make an adiabatic assumption considering small signal 
frequencies that makes the NEP valid at each time for the corresponding value of the signal. 

Using Eq. (|2l|) we can integrate Eq. (^) with the initial condition xit^) = xo and obtain 
the conditional probability ni{t \ xq, to)- This result allows us to calculate the autocorrelation 
function, the power spectrum and finally the SNR. The details of the calculation were shown 
in Ref. 0,|16[- When the symmetrical case is considered all the results reduce to those in 



T^ . For the SNR, up to the relevant (second) order in the signal amplitude A, we find for 



the SNR the result 16 



4/ii/i2(/Ui +/i2) 
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In order to evaluate the transition rates between both states we discretize the space and 
the fields as 



X Xi, v{x)) (Ui, U2..., UN, Vi, 



(23) 



and use the Kramers-like formula 15 



W, 



exp 



(24) 



2vr^ \<I>'L\ ' L V 
where A+ is the unstable eigenvalue of the deterministic fiux at the unstable state Um, 
and $^ indicate the determinants of the matrix of second order derivatives of the NEP with 
respect to the discretized fields in the states Ui and Um respectively, and and are the 
values of the NEP evaluated at the stationary states Ui and Um, i = 1,2. Finally, in order 
to compute the SNR as indicated above, we calculate the parameters /i, and ctj numerically 
as 

dWi 



f^i = Wi\s{t). 



■0 ) 



on 



dS{t) 



(25) 



5(i)=0 



B. Role of the NEP symmetry 

To fix ideas we consider the system given by Eq. (|16|) in one dimension, with the spatial 
coordinate x varying between —L and L, and assuming Dirichlet boundary conditions for 
the three fields. We adopt the following piecewise linealized form for the nonlinear function 

fin) 

f{u) = -u + 9{u - a) (26) 

where 9{u) is the step function. We fix = 0, g2 = .02, = .01 and g2 = 1. This leads us 
to a situation in which we have essentially only one noise source (^2) acting on v. With this 
choice we have Qu <^ Qv and 7 results approximately independent of /?. 

The piecewise linearized form of f{u) allows us to calculate analytically the stationary 
inhomogeneous patterns of the associated deterministic system as linear combinations of 
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exponentials |]8|JT0[]. We search for solutions symmetric with respect to x = 0, with a central 
activated region {u > a) surrounded by a non activated region {u < a). Depending on the 
values of a, D, (3 and (3' we find four (two stable and two unstable), two (one stable and one 
unstable) or zero stationary inhomogeneous solutions. Furthermore the homogeneous null 
solution {u = V = w = always exists and is stable. We call these solutions Si and nsi 
(existing in the regions were there are two or four inhomogeneous solutions), S2 and ns2 
(existing only in the region of four solutions) and sq (the null solution). The Sj's are the 
stable solutions and the nsj's are the unstable ones. 

We will focus our analysis on a region of parameters where the system has two stationary 
stable patterns (stationary lineary stable solutions of Eqs. (|12D for -u, f, and w) and one 
stationary unstable pattern (stationary lineary unstable solution of Eqs. (0)). In Fig. 3 



of Ref. [|T2| the w-fields for the three patterns for a particular choice of the parameters were 
shown. We call Ui the large stable pattern which has a central activated region {u > a), 
U2 the small stable pattern which reduces to the homogeneous null solution when S is set 
equal to zero, and Um the unstable pattern. A complete study of the pattern formation of 
this system can be found in Ref. |T2I . 



In the region of only two stable patterns we are considering, the deterministic dynamics 
given by Eqs. ([16|) drives the system toward one of the patterns (selected depending on the 
initial condition) which is reached asymptotically. If small fiuctuations are present in the 
system the fields fiuctuate around one of the stable patterns and transitions between the 



two patterns become possible. Note that the (yff in Eqs. (|T6D are constants that couple the 
noises to the system while the intensity of the fiuctuations is determined by the parameter 
V- 

It is worth mentioning that the NEP given in Eq. (|18|) is valid for the system in Eq. 
(p!6|) for arbitrary number of spatial dimensions, for an arbitrary nonlinear function f{u), 
and with the parameter region of validity being independent of the choice of f{u) ||I2|. The 
consideration of only one spatial dimension and the particular election of f{u) are in order to 
simplify the calculations, particularly regarding pattern formation. The signal is introduced 
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as a (slow) modulation in a parameter 5* by setting S = S(t) = A cos^Ugt) added to the 
autocatalytic function f{u). 

We now analyze the SR phenomenon in our spatially extended system using the theory 
discussed before. To proceed with such an analysis we identify the two stable patterns (f/i 
and U2) with the states of the two-state-theory. Hence, the discrete variable x will adopt 
values Ci and C2 according to whether the system is in the states Ui or U2, yielding the result 
for the SNR in Eq. (p2|). The changes induced in the patterns and their stability by the 
variation of some model parameter will be reflected in changes in the values of /ij and 
and, accordingly, will affect the results for the SNR. 

We fix L = = /5' = 1,7 = 10.026,7' = z/ = 10,< = l,g^ = 0,g^ = .05 and 
(^2 = -01, and leave D, a and rj (the noise intensity) as free parameters. Note that with the 
chosen values for the (7f' s, the only relevant noise term in the system (Eq. (0)) is 5'i^i(t) in 
the equation for u that appears added to the signal (hence it can be considered as coming 
together with the signal). The parameters and g2 are set different from zero to keep 
the system inside the parameter region where as defined in Eq. ( [T^ ) is valid as a 

nonequilibrium potential [|I^]. 

In Fig. 2 we show the results for the SNR (R) as a function of the noise intensity for 
different values of D and a. We see that while keeping a constant (a = .25) (Fig. 2 a), the 
largest values of R are those for D = Dg, which is the symmetric situation. Also, if we fix 
D = Dg (Fig. 2 b), any departure of a from the value .25 (that is any departure from the 
symmetric situation) lowers the values of R. Hence, the symmetric situation is found to be 
the most favorable one concerning the improvement of SNR. Note that the maximum of the 
R vs. rj curve (for fixed values of a and D), that we will call Rmax, increases with symmetry 
and reaches its largest value for the symmetric situation. In Fig. 3 we show Rmax plotted 
as a function of D for a = .25, where it is apparent that the optimum value of diffusion is 
D = Dg corresponding, as indicated, to the symmetric case. 

A fact that arises from these results is that, while keeping all the other parameters of the 
system fixed, there exists an optimal value of diffusion (coupling of the distributed system) 
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that maximizes SNR. The interesting aspect is that this optimal value is the one that makes 
the potential symmetric. 

It is worth mentioning that these results do not contradict but complete those in ||^ where 
enhancement due to coupling was found, since in that work only symmetric situations were 
analyzed. Roughly speaking, the main result in can be summarized saying that, given 
two different symmetric situations (each one necessarily having different values of D and 
a), the one with the higher value of D produces higher values of SNR. However, we must 
keep in mind that for a too large value of some of the approximations involved in the 
calculations may break down |^. Also, it is worth here pointing out that the same thing 
may happen for too large asymmetries. For example, consider an extremely asymmetric 
situation where the barrier for, say the transition from state Ui to t/25 is much larger than 
the barrier for the opposite transition. In such a case, the values of the noise intensity 
leading to reasonable jumping rates from Ui to U2 will be far beyond the validity of the 
Kramers-like approximation for the inverse transition. 

C. Enhancement due to selective coupling 

The model we discuss next is an extension of the one-component RD model discussed 
in Section 2, but now the diffusive parameter depends on the field (f){x,t). As a matter of 
fact, since in the ballast resistor the thermal conductivity is a function of the energy density, 
the resulting equation for the temperature field includes a temperature-dependent diffusion 
coefficient in a natural way. 

By adequate rescaling of the field, space-time variables and parameters, we get a dimen- 
sionless time-evolution equation for the field t) 

dt4>{x,t) = {D{^)d,^) + /(</.) + C{x,t), (27) 

where ^(x, t) is a white noise in space and time, and /(</>) = —cj) + 6{(j) — (pc), is the step 
function. All the effects of the parameters that keep the system away of equilibrium (such 
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as the electric current in the electrothermal devices or some external reactant concentration 
in chemical models) are included in 0^. 

As was done for the reaction term [^,^|10|, a simple choice that retains however the 



qualitative features of the system is to consider the following dependence of the diffusion 
term on the field variable 

= D,{l + hd[ct>-ct>,)l (28) 

For simplicity, here we choose the same threshold 0c for the reaction term and the diffusion 
coefficient. 

We assume the system to be limited to a bounded domain x G [—L^ L] with Dirichlet 
boundary conditions at both ends, i.e. (j){±L,t) = 0. As before, the piecewise-linear approx- 
imation of the reaction term in Eg. ([27|) was chosen in order to find analytical expressions 
for its stationary spatially-symmetric solutions. In addition to the trivial solution (poi^) = 
(which is linearly stable and exists for the whole range of parameters) we find another lin- 
early stable nonhomogeneous structure (psi^) — presenting an excited central zone (where 
0s (x) > 0c) for —Xc < X < Xc — and a similar unstable structure 0u(x), which exhibits a 
smaller excited central zone. The form of these patterns is analogous to what has been 



obtained in previous related works ||8|,[10|1. The difference is that in the present case d(l)/dx\x^ 
is discontinuous and the area of the "activated" central zone depends on h. 

The indicated patterns are extrema of the NEP. In fact, the unstable pattern 0u(x) is a 
saddle-point of this functional, separating the attractors (po{x) and 0s (x). For the case of a 
field-dependent diffusion coefficient D{(j){x,t)) as described by Eq. (|27|) , the NEP reads 

$[0] = |_^^^|-|^%(0')/(0')rf0' + ^/^(0)|^) jrfx. (29) 

Given that dt4> = — (l/-D(0))5$/50 one finds $ = — / [5^ / dcj))^ dx < 0, thus warranting 
the LF property. 

For a given threshold value 0*, both wells corresponding in a representation of the NEP 
to the linearly stable states have the same depth (i.e. both states are equally stable). Figure 
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4 shows the dependence of $[0] on the parameter 0c- As in previous cases, we analyze only 
the neighborhood of 0c = 0* HI- Here we also consider the neighborhood oi h = 0, where 
the main trends of the effect can be captured. 

In Fig. 5 we depict the dependence of R on the noise intensity 7, for several (positive) 
values of h. These curves show the typical maximum that has become the fingerprint of the 
stochastic resonance phenomenon. Figure 6 is a plot of the value Rmax of these maxima as 
a function of h. The dramatic increase of Rmax (several dB for a small positive variation of 
h), is apparent and shows the strong effect that the selective coupling (or field- dependent 
diffusivity) has on the response of the system. 

It must be noted that the only two approximations made in order to derive 
Eq.(^2]) — namely the Kramers-like expression in Eq.(p^) and the two-level approximation 
used for the evaluation of the correlation function 0Jl^ — break down for large positive val- 
ues of h because for increasing selectivity the curves of $[0] vs. 0c in Fig. 4 shift towards the 
left, which in turn means that the barrier separating the attractors at 0* tends to zero. This 
effect is basically the same as the one discussed in Ref. |^ in connection with global diffu- 
sivity Dq. It is also worth noting that except for the two aforementioned approximations, 
all the previous results (e.g. the profiles of the stationary patterns and the corresponding 
values of the non-equilibrium potential) are analytically exact. 



IV. CONCLUSIONS 

In this work we have presented the basic ingredients for studying SR in coupled or ex- 
tended media {stochastic resonant medium |^) by means of nonequilibrium potential tech- 
niques. We have shown some examples of NEP for RD systems and the way we can exploit 
them to get information about the system's response (the SNR) to weak potential modula- 
tions. 

In previous works we have shown, in agreement with numerical simulations |^ how the 
coupling enhances the system's response. Here we have shown two other aspects: (a) the 
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role played by the symmetry of the NEP on such enhancement and; (b) the effect of selective 
coupling. 

In the first case the results we have obtained clearly indicate the central role played by 
symmetry in improving the SNR. We studied the behavior of Rmax, that is the maximum 
of the SNR vs. rj curve, as the different model parameters are (not simultaneously) var- 
ied, finding that Rmax always increases with the symmetry of the potential. This fact led 
us to our main result: the optimal values of the different model parameters (for instance 
diffusivity or threshold), as regards the maximization of Rmax, correspond to those making 
the potential more symmetric in each situation. Besides the analysis of the influence of 
symmetry on stochastic resonance, it is important to remark that the mere consideration of 
asymmetric situations has its own relevance. This is because such bistable asymmetric mod- 
els provide, for example, the appropriate framework for describing SR in voltage-dependent 
ion channels, as proposed in some biological experiments. In those systems, the conducting 
state is associated a higher-energy well than the non-conducting one [|T^. 

In the second case, the results of our analysis of a scalar RD system (a generalization of 
the continuous limit of the model analyzed in 0,^) indicate that a "selective" coupling (that 
is a field diffusion dependent constant) could offer a new enhancing mechanism in coupling 
systems. This prediction prompts to devise experiments (for instance, through electronic 
setups) as well as numerical simulations taking into account the indicated selective coupling. 

One direction in which the present studies can be extended corresponds to the analysis 
of the competition between local and nonlocal couplings. Consider the same NEP described 



by Eq. (|T8D, but with only one inhibitor, the temporally slaved w. For this system the NEP 



$ has the form 

t), v{x, t)] = Jdx [— (Vm(x, t) f + C J dx'G{x , x')u{x,t)u{x' ,t)\ , (30) 

where the coupling constant C, that in one of the systems studied in was positive (that is 
due to an inhibition interaction), is now assumed negative (that is an excitatory interaction) . 
The results of preliminary studies for the case of excitatory coupling, where we used the 
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same kind of approach indicated above, are shown in Fig. 6. There we have compared the 
results for the SNR vs. noise intensity for the FitzHugh-Naguno model |^ with inhibitory 
interaction, with those obtained in the above indicated case. It is apparent, on one hand, the 
SNR enhancement due to the excitatory interaction, while on the other, that small changes 
in the diffusion constant of the activator ("local coupling" constant) strongly affect the SNR. 
The analysis of the different aspects of this problem will be the subject of a forthcoming 
paper Ijl^ . 



Acknowledgments: The authors want to thank V. Griinfeld for a revision of the 
manuscript. HSW thanks to F.Castelpoggi, R.Deza, G.Drazer, G.Iziis, M.Kuperman, 
D.Zanette, for their collaboration in early stages of this work; to R. Graham, E. Tirapegui 
and R. Toral for useful discussions; and acknowledges financial support from CONICET 
(Project PIP-4953/96). 



18 



REFERENCES 



[1] M.C. Cross and P.C. Hohenberg, Rev. Mod. Phys. 65, 851 (1993); E. Meron, 
Phys. Rep. 218, 1 (1992), J.D. Murray Mathematical Biology (Springer- Verlag, 1985); 
A.S. Mikhailov, Foundations of Synergetics I, (Springer- Verlag, 1990); D. Walgraef, 
Spatio-temporal pattern formation, (Springer- Verlag, New York , 1997); 

[2] W. Horsthemke and R. Lefever, Noise induced transitions, (Springer- Verlag, New York, 
1986); C. Van den Broeck, J. M. R. Parrondo, R. Toral and R. Kawai, Phys. Rev. 
E 55, 4084 (1995); S. Mangioni, R. Deza, R. Toral and H.S. Wio, Phys. Rev. E 61, 
223 (2000); P. Hanggi and R. Bartussek, Brownian Rectifiers in: Nonlinear Physics of 
Complex Systems, J. Parisi, S. C. Miiller, and W. Zimmermann (Eds.), Lecture Notes 
in Physics 476, 294 - 308 (1996); for an update review on Brownian motors see P.T. 
Reimann, Brownian motors: noise transport far from equilibrium, to appear in Physics 
Reports; P.T. Reimann, R. Kawai, C. Van den Broeck and P. Hanggi, Europhys. Lett. 
45, 545 (1999); S. Mangioni, R. Deza and H.S. Wio, Phys. Rev. E 63, 041115 (2001); 
J. Buceta, J.M. Parrondo, C. Van den Broeck and F.J. de la Rubia, Phys. Rev. E 61, 
6287 (2000); R. J. Deissler, J. Stat. Phys. 54, 1459 (1989); K. L. Babcock, G. Ahlers 
and D. S. Cannell, Phys. Rev. E 50, 3670 (1994); M. Neufeld, D. Walgraef and M. San 
Miguel, Phys. Rev. E 54, 6344 (1996). 

[3] L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, Rev. Mod. Phys., 70, 223 
(1998). 

[4] J. Lindner, B. Meadows, W. Ditto, M. Inchiosa and A. Bulsara, Phys. Rev. Lett. 75, 3 
(1995); ibid, Phys. Rev. E 53, 2081 (1996). 

[5] H. S. Wio, Phys. Rev. E 54, R3075 (1996); H.S. Wio and F. Castelpoggi, Proc. Conf. 
UPoN'96, C.R.Doering, L.B.Kiss and M.Schlesinger Eds. (World Scientific, Singapore, 
1997), pg. 229; F. Castelpoggi and H.S. Wio, Europhysics Letters 38, 91 (1997); F. 
Castelpoggi and H. S. Wio, Phys. Rev. E 57, 5112 (1998); M. Kuperman, H.S. Wio, G. 

19 



Iziis and R. Deza, Phys. Rev. E 57, 5122 (1998); M. Kuperman, H.S. Wio, G. Izus, R. 
Deza and F. Castelpoggi, Physica A 257, 275 (1998). 

[6] R.Graham, Weak Noise Limit and Nonequilibrium Potentials of Dissipative Dynamical 
Systems, in Instabilities and Nonequilibrium Structures, Eds.E.Tirapegui and D.Villaroel 
(D.Reidel, Dordrecht, 1987); R. Graham and T. Tel, Phys. Rev. A 42, 4661 (1990); 
O.Descalzi and R.Graham, Phys. Lett. A 170, 84 (1992); O.Descalzi and R.Graham, 
Z.Phys.B 93, 509 (1994). 

[7] R.Montagne, E.Hernandez-Garcia and M.San Miguel, Physica D 96, 47 (1996). 

[8] G. Izus, O. Ramirez, R. Deza, H.S.Wio, D.Zanette and C.Borzi, Phys.Rev.E 52, 129 

(1995) ; D.H.Zanette, H. S. Wio, and R. Deza, Phys. Rev. E 53 353 (1996); G. Izus, 
H. S. Wio, J. Reyes de Rueda, O. Ramirez and R. Deza, Int. J. Mod. Phys. B 10, 1273 

(1996) ; G. Izus, R. Deza, C.Borzi, and H.S.Wio, Physica A 237, 135 (1997); G. Izus, R. 
Deza, H.S.Wio, and C.Borzi, Phys.Rev.E 55, 4005 (1997); P. Castelpoggi, H.S.Wio, and 
D.H. Zanette, Int.J. Mod.Phys. B 11, 1717 (1997); G. Drazer and H.S. Wio, Physica A 
240, 571 (1997). 

[9] H.S. Wio, Nonequilibrium Potential in Reaction- Diffusion Systems, in: Fourth Granada 
Lectures in Computational Physics, editores P.L. Garrido and J. Marro, (Springer- 
Verlag, 1997), pg.l35. 

[10] G. Izus, R. Deza, O. Ramirez, H. Wio, D. Zanette and C. Borzi, Phys. Rev. E 
52, 129 (1995); G. Iziis, H. S. Wio, J. Reyes de Rueda, O. Ramirez and R. Deza, 
Int.J.Mod.Phys. B 10, 1273 (1996). 

[11] D. M. Petrich and R. E. Goldstein; Phys. Rev. Lett. 72, 1120 (1994); R. E. Goldstein, 
D.J. Muraki and D. M. Petrich; Phys. Rev. E. 55, 3993 (1996). 

[12] S. Bouzat and H.S. Wio, Phys. Lett. A. 247, 297 (1998). 

[13] G. Iziis, R. Deza, and H.S.Wio, Phys. Rev. E 58, 93 (1998); G.G. Izus, R.R. Deza and 



20 



H.S. Wio, Comp. Phys. Comm. 121-122, 406 (1999). 

[14] C. Nicolis, Tellus 34, 1 (1982); B. McNamara, K. Wiesenfeld, Phys. Rev. A 39, 4854 
(1989). 

[15] P. Hanggi , P. Talkner and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990) 

[16] F. Castelpoggi and H.S. Wio, Physical Review E 57, 5112 (1998); H.S.Wio and S. 
Bouzat, Brazil. J. Phys 29, 136 (1999); S. Bouzat and H.S. Wio, Phys. Rev. E 59, 5142 
(1999); B. von Haeften, R. Deza and H. S. Wio, Physical Review Letters 84, 404 (2000). 

[17] S. M. Bezrukov and I. Vodyanoy, Nature 378, 362 (1995); I. Goychuk and P. Hanggi, 
Phys. Rev. E 61, 4272 (2000); A.D. Sanchez, J.A. ReveUi and H.S. Wio; Phys. Lett. A 
277, 304 (2000). 

[18] B. von Haeften and H. S. Wio, in preparation. 



21 



FIGURES 

FIG. 1. NEP evaluated at the stationary patterns, for Dirichlet b.c, as a function of for 

L = 1 and D = 1. The bottom curve (1) corresponds to (ps{x) and the top one (3) to 4>u{x). The 
bistabiUty point is indicated. 

FIG. 2. a) SNR as a function of the noise intensity for a = .25 and different values of D. The 
solid line corresponds to the symmetric situation D = Dg, the long-dashed line to D = .35 and 
the short-dashed line to D = .25. b) SNR as a function of the noise intensity for D = Dg and 
different values of a. The solid line corresponds to the symmetric situation a = .25, the dotted line 
corresponds to a = .27 and the dot-dashed to a = .23. 

FIG. 3. Maximum of SNR (Rmax) as a function of the activator diffusion D for a = .25. The 
maximum of Rmax occurs for the symmetric situation D = Dg 

FIG. 4. SNR (R) as a function of the noise intensity (here indicated by 7), for three values of 
h: h = 0.0 (fuh line), h = -0.25 (dashed line) and h = 0.25 (dotted line). We have fixed L = 1, 
DQ = l,6(t)c = 0.01 and n = 0.01. 

FIG. 5. Maximum R^ax of the SNR curve (Fig.4) as a function of h, for three values of Dq: 
Dq = 0.9 (dashed line), Dq = 1. (full line) and Dq = 1.1 (dotted line). The arrows a and b 
indicate the response gain due to an homogeneous increase of the coupling and to a selective one 
respectively. The larger gain in the second case is apparent. The inset shows the dependence of 
Rmax on Do for h = —0.25 (lower line), h = and h = 0.25 (upper line). 

FIG. 6. SNR for the case of competition between local and nonlocal interaction as a function 
of the noise intensity (here indicated by 7). All curves have Dy = 4 (nonlocal interaction). The 
following curves correspond to a nonlocal excitatory interaction with different values of the 
upper broken line Du = 1.2, the continuous line Du = 0.8, the lower broken line Du = 1. The 
dotted line corresponds to the original (inhibitory) FitzHugh-Nagumo case = 1.2). 
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